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Vortex core reversal in magnetic particle is essentially influenced by a surface anisotropy. Under 
the action of a perpendicular static magnetic field the vortex core undergoes a shape deformationof 
pillow- or barrel-shaped type, depending on the type of the surface anisotropy. This deformation 
plays a key point in the switching mechanism: We predict that the vortex polarity switching is 
accompanied (i) by a linear singularity in case of Heisenberg magnet with bulk anisotropy only and 
(ii) by a point singularities in case of surface anisotropy or exchange anisotropy. We study in details 
the switching process using spin-lattice simulations and propose a simple analytical description 
using a wired core model, which provides an adequate description of the Bloch point statics, its 
dynamics and the Bloch point mediated switching process. Our analytical predictions are confirmed 
by spin-lattice simulations for Heisenberg magnet and micromagnetic simulations for nanomagnet 
with account of a dipolar interaction. 
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I. INTRODUCTION 

Magnetization reversal in small magnetic particles 
is one of the fundamental issues of the modern mag¬ 
netism. Different concepts of switching, including switch¬ 
ing by magnetic fields (quasistatic and precessional ones), 
switching by spin injection and all optical switching 
are intensively studied and widely used for applications 
in magnetic data storage. Nowadays inhomogeneous 
switching,^ i. e. the process of magnetization rever¬ 
sal for inhomogeneous magnetization configurations, at¬ 
tracts growing interest. Particular attention is paid to 
topologically protected configuration like magnetic vor¬ 
tex. The vortex is determined by a planar closed flux- 
free configuration in a sample with a localized nonplanar 
magnetization, so called vortex core. The sense of the 
core magnetization direction is characterized by the vor¬ 
tex polarity (up, p = -|-1 or down, p = —1). Because of 
the topological stability of vortices with different polari¬ 
ties, one can consider the vortex polarity as a bit of in¬ 
formation in nonvolatile magnetic vortex random-access 
memories.The corresponding inhomogeneous magne¬ 
tization reversal in the vortex-state particles is called the 
vortex core reversal or the vortex polarity switching. 

There exist several concepts of the vortex polarity 
switching, e. g. by magnetic fields, by spin injec¬ 
tion, see Ref. 7 and references therein. Independently 
of these concepts one can separate two main scenar¬ 
ios: axially-symmetric (or punch-through) scenario and 
axially-asymmetric one. In the first scenario the vor¬ 
tex does not make the macroscopic motion, the vortex 
core reversal is caused by the direct excitation of radially 
symmetric magnon modes. Such a switching occurs, e. g. 
under the action of dc transversal magnetic field®^^^, or 
by ac transversal fields^^^^®. The axially-asymmetric sce¬ 
nario is caused by nonlinear resonance between coupled 
magnon modes;^^’^® it is accompanied by the temporary 
creation of the vortex-antivortex pairs^®’^°. Such a pic¬ 


ture was observed in different setups.' 

Here we consider the axially-symmetrical switching of 
the vortex polarity. During the reversal process the main 
changes occur in a circular Bloch line^^, which is the line 
passing trough the sample thickness and connecting cen¬ 
ters of all vortex cross-sections. The magnetization direc¬ 
tion in the Bloch line corresponds to the vortex polarity: 
it has to reverse its direction during the switching pro¬ 
cess. In case when a reversal occurs simultaneously for 
the whole Bloch line 7 , there appears a line singularity 
during switching: the magnetization vector m vanishes 
at the switching moment for the whole line, m{'y) = 0 , 
and the temporary pure planar vortex (linear singularity) 
appears during the switching. Such a two-dimensional 
(2D) switching process was considered in Ref. 11 to de¬ 
scribe the reversal in very thin magnets. Another possi¬ 
bility was suggested by Arrott et al. . Arrott’s model is 
stray-field-free, where the switching is made possible by 
propagating of point singularities, Bloch points with co¬ 
ordinates In that case the reversal occurs non- 

simultaneously: the magnetization vanishes only in one 
or a few points, the so-called Bloch points, m(rBp) = 0, 
and the switching is accompanied by the nucleation, mo¬ 
tion and final annihilation of the Bloch points. Such a 
process was studied by Thiaville et al. using micromag¬ 
netic simulations: the switching process is typically ac¬ 
companied by the creation of two Bloch points, however, 
the single Bloch point scenario was also mentioned. 

The purpose of the current study is to describe process 
of the axially-symmetrical switching, to find the mecha¬ 
nism, which is responsible for simultaneous or nonsimul- 
taneous switching, to understand how many Bloch points 
are nucleated during the switching, and to describe this 
process analytically. 

Point singularities were introduced in magnetism by 
Feldtkeller and Doring^®. For example, in the vicinity 
of the Bloch point, situated at rpp = 0, the magne¬ 
tization distribution has a hedgehog-like configuration. 
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m = ±Qrlr with Q being the constant orthogonal 
matrix.^® They are analogues of magnetic monopoles in 
elementary particle physics for spin waves.Bloch points 
were observed as twisting of the Bloch lines in garnet 
crystals^®’^® and as singular point of the Bloch point do¬ 
main wall in magnetic nanowires^®. According to sim¬ 
ulations the Bloch points appear as transient states in 
the vortex dynamics under the action of dc^° and ac^^^s 
transversal fields; as temporal formations (micromag- 
netic drops) during the magnetization reversal in soft 
magnetic cylinders.The Bloch points are also created 
during the vortex-antivortex annihilation processes. 

Theoretical treatment of the Bloch point is a compli¬ 
cated task. Even a “simple” problem of the Bloch point 
structure in a spherical particle caused a long-time dis¬ 
cussion, see Ref. 33 and references therein. The problem 
becomes complicated due to the fact that the Bloch point 
distribution with m(rBp) = 0 does not fulfil the condi¬ 
tion of magnetization normalization, \m\ = 1. There 
is no physical singularity in reality. Due to the dis¬ 
creetness of the spin lattice, the Bloch point (as well as 
linear singularity) is situated in an interstitial site po¬ 
sition, hence the singularity appears only in continuum 
description. That is why the modeling of Bloch points is 
also nontrivial task: there appear mesh dependent prob¬ 
lems, e.g. mesh-friction effect and a strong mesh depen¬ 
dence of the switching field.Therefore specific simu¬ 
lators with atomistic resolution are needed to overcome 
such difficulties.The role of discreetness of the mag¬ 
netic lattice becomes of great importance; in particular, 
namely the lattice creates an effective pinning field for 
the Bloch point;^^ it changes the Bloch point behavior 
relatively to the main crystallographic directions when a 
depinning held for one direction can be much lower than 
in other direction. 

The problem of dynamics of Bloch points is not less 
challenging. The spectrum of Bloch point oscillations 
along the Bloch line was measured for yttrium garnet 
ferrite.^® Dynamical properties of the Bloch point, its 
mobility and mass were derived for the Bloch point in¬ 
side the Bloch line.^^^"^° We are interesting in a specihc 
dynamics the Bloch point, namely, the Bloch point medi¬ 
ated switching of the vortex polarity. It is already known 
from our recent studies^^^s dominating contri¬ 

bution to the switching mechanism is caused by the ex¬ 
change interaction inside the vortex core. Therefore one 
can expect that the dipolar interaction, which is essen¬ 
tial for the statics and slow dynamics of the magnet, in 
particular, for the vortex state nanodot, does not play 
the dominant role for the switching process. 

Argued by above mentioned reasons of the dominant 
role of exchange interaction we do not take into ac¬ 
count in most cases the dipolar interaction. Instead 
we additionally include into consideration the surface 
anisotropy effects. In a disk-shaped particles the surface 
anisotropy favors magnetization curling states of onion 
(two half vortex^^ or capacitor^^) configuration and a 
vortex one^^’^^. Below we consider the vortex configura¬ 


tion. One has to note that the vortex is the metastable 
state in model under consideration,^^ while it can form 
the ground state in models with finite skin depth of the 
surface anisotropy'*^^ or in models with dipolar induced ef¬ 
fective inhomogeneous anisotropy.In the last case there 
appear two effective inhomogeneous anisotropy terms: 
one is effective easy-plane anisotropy of face surface 
charges and another one is effective anisotropy of edge 
surface charges;^^ namely these two effective anisotropies 
are responsible for the the vortex state conhguration. 

In the current study we consider the Heisenberg disk¬ 
shaped magnet in a vortex configuration. Without sur¬ 
face anisotropy the vortex prohle does not depend on a 
thickness coordinate z, it has a 2D shape. The pres¬ 
ence of the surface anisotropy breaks such a symme¬ 
try and results in deformation of the vortex core pro¬ 
file. Depending on the type of surface anisotropy, the 
vortex profile becomes barrel- or pillow-shaped deformed 
for the easy-surface (ES) and easy-normal (EN) surface 
anisotropies, respectively.^^ We will see below that in¬ 
homogeneity of the Bloch line is crucial for understand¬ 
ing the Bloch point mediated switching mechanism: the 
switching process in Heisenberg three dimensional (3D) 
magnet with bulk anisotropy only is mediated by the 
linear singularity. In contrast, the switching in magnet 
with surface anisotropy is accompanied by the simulta¬ 
neous nucleation of two Bloch points on the face surfaces 
of the sample for the EN case and by the nucleation of 
two Bloch points inside the sample (at the center of the 
Bloch line) for the ES case. We also analyze the case 
of the exchange anisotropy; qualitatively the mechanism 
is very similar to the EN surface anisotropy: during the 
switching two Bloch points enter the sample on face sur¬ 
faces, move inside to each other and finally annihilate at 
the center of the disk. We study in details the process of 
the vortex core reversal, compute the Bloch points tra¬ 
jectories, and hnd how the Bloch point speed depends 
on the sample thickness. We propose a simple analytical 
picture which describes both statics of the vortex in the 
cone phase with account of the surface anisotropy effects 
and dynamics of the vortex polarity switching process, 
which is in agreement with the full-scale spin-lattice sim¬ 
ulations. We also performed micromagnetic OOMMF^® 
simulations with account of surface anisotropy and ad¬ 
ditional dipolar interaction, which confirm main features 
of our vortex core deformation. 

The paper is organized as follows. In the Sec. H we de¬ 
scribe the model of the Heisenberg magnet with account 
of both bulk and surface anisotropies. In the Sec. HI 
we present results of spin-lattice simulations of vortex 
polarity switching for Heisenberg magnets with bulk 
easy-plane single-ion anistropy (Sec. HI A), with addi¬ 
tional easy-surface and easy-normal surface anisotropies 
(Sec. HIB) and with bulk exchange anisotropy (Ap¬ 
pendix B). We make analytical description of the switch¬ 
ing phenomenon in the framework of proposed wired core 
model which demostrates a good qualitative agreement 
with simulations in Sec. IV. To validate our results for 
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nanoparticles we modeled a vortex state nanodisk us¬ 
ing OOMMF micromagnetic simulations with account of 
dipolar interaction and surface anisotropy, which are dis¬ 
cussed in Sec. V. In Appendix A we consider how the 
homogeneous state is deformed with account of the sur¬ 
face anisotropy. The basic equations for the wired core 
model are derived in Appendix D. 


Here t = a\f\ is the magnetic length, a is the lattice con¬ 
stant and A = J/K, xr = Ks/K is the surface anisotropy 
measured in units of K and is the normal to the sur¬ 
face. We will consider cases of small surface anisotropy 
\x\ < 1. The parameter h is the field intensity, normal¬ 
ized by the anisotropy held Ha, 

h = H/Ha, Ha = 


II. THE MODEL AND THE VORTEX 
SOLUTION 

Let us consider a classical lattice Heisenberg magnet 
with a simple cubic lattice and Hamiltonian 


The equilibrium magnetization distribution can be found 
by variation of the energy functional (3) which results in 
the following boundary-value problem:^^’^^ 

m X + h — {m ■ z)z'j = 0, (4a) 


df = - J§^ ^ rrin ■ m-n+s - • z) + J{“, 

(ti. 5) n 

(la) 

where J > 0 is the exchange integral, S = 
1/2 is the length of a classical spin, Tn„ = 
{^/T^^m^cos (pn, is the normal¬ 

ized magnetic moment on a 3D site position n, the 3D 
index S runs over the nearest neighbors, is the Bohr 
magneton, H is the intensity of external dc magnetic 
held, directed along z-axis, and is the Hamilto¬ 
nian of single-ion anisotropy. The latter includes the 
bulk easy-plane anisotropy (with the anisotropy constant 
AT > 0), which favors magnetization distribution within 
xy-plane, and the Neel surface anisotropy term (with the 
anisotropy constant KgY^, 

^ ^(m„ • zf - {mi ■ uisf. (lb) 

n (l,S) 


where index I runs over the surface sites and the unit 
vector Uis connects the nearest neighbors of the lattice. 
Dynamics of magnetization is governed by the discrete 
version of the Landau-Lifshitz-Gilbert equation 


dmn dJ^f dm„ 

—— = m„ X -h emn x ——, 

dr omn dr 


( 2 ) 


where r = Dt is the normalized time, D = KB/h, 
Jff = df/(ATS^) is the normalized energy, h is the Plank’s 
constant, and e is the Gilbert damping constant. 

The continuum approach is based on smoothing of the 
lattice model using the normalized magnetization 


m(r, t) = (^\/l 


— m? cos (/}, \/1 — sin < 





K 


(m • rig) [(m • ns)m — ris] 


S 


(4b) 


Without the surface anisotropy {x = 0) the ground state 
of the uniaxial magnet is determined by the field in¬ 
tensity. As it follows from (4a), the magnetization z- 
component m = h = const, the boundary condition (4b) 
has the Neumann type, dm/dria = 0 which is trivially 
satisfied for the constant m solution. 

The presence of the surface anisotropy drastically 
changes the symmetry of the problem. The boundary 
conditions (4b) becomes of the Robin type"*’®, which is 
the source of the symmetry breaking of the magnetization 
structure and causes the effect of magnetization curling. 
There appears nonhomogeneous ground states with typ¬ 
ical examples as hedgehog, throttled, artichoke, onion, 
vortex configurations, 

Let us start with the easy-plane magnet which has the 
shape of infinite him of a hnite thickness L = {N — l)a, 
where N is the number of lattice sites along thickness 
coordinate z. Without external held and the surface 
anisotropy the ground state of such a magnet is an uni¬ 
form easy plane magnetization, m = 0. The same state 
takes place with account of the weak surface anisotropy 
in the no-magnetic held case. 

The magnetization conhguration under the action of 
applied held is affected by the surface anisotropy inhu- 
ence: the magnetization distribution breaks its homo¬ 
geneity in z-direction due to the surface anisotropy inhu- 
ence 


mo(z, h) = h 


1 


— X' 


1 — cosh(-\/l — h?z/C) 
A sinh(-\/l — h‘^L/2£) 


( 5 ) 


where m = m{r,T) and p = 4 >{r,T). The total energy, 
which corresponds to the Hamiltonian (1), takes the form 


a :§2 

A=^(4 + ^s) 

= / dU 


-m • V m H- m ■ h 

2 2 


( 3 ) 


S’s = ^ J dS{m ■ rig)'^. 


see Appendix A for details. Typical magnetization prohle 
mo as function of z for EN and ES surface anisotropies is 
shown in Fig. 1. Symbols correspond to the spin-lattice 
simulations,®^ see Sec. HI for details. One can see that 
the sign of the curvature mo = f{z) is determined by the 
sign of the product xh. 

In the following paper we consider the disk¬ 
shaped samples with such parameters, that the surface 
anisotropy favors the vortex conhguration. ^**’^2 
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ES, sim. • ES, an. — EN, sim. ■ EN, an. — 



FIG. 1: (Color online) Magnetization profile mo(z;/i) as 
a function the thickness coordinate z foi h = 0.7 (level 
is shown by dotted line) and |/c| = 0.25. Symbols 
correspond to spin-lattice simulations and solid lines to 
the first order expansion of Eq. (5). Schematics in 
insets show the magnetization distribution in the film 
with surface anisotropy for the both cases (red arrows 
show magnetization direction and blue arrows show 
field direction). 


plane structure of the stationary vortex state (j){r) is de¬ 
termined by the following magnetization distribution 

</> = X + C'|, (6a) 

where (r, %, z) are the cylinder coordinates. Without sur¬ 
face anisotropy effects, the vortex chirality C can take 
any value due to the isotropy of Heisenberg exchange, 
while the surface anisotropy on the edge surface fixes its 
value: C = ±1 for ES magnets (>f > 0) and C = 0 or 
C = 2 for EN magnets (xr < 0).^'*’ 

Let us discuss the vortex out-of-plane structure, m(r). 
Without magnetic field the out-of-plane vortex profile has 
a bell-shaped structure^^ 

m(r,z)«pexp(^-^^^), (6b) 

which generalizes a well-known Feldtkeller Ansatz,^^’®^ 
originally used for magnets without surface anisotropy, 
when the dimensionless core width ic « 1. Depending 
on the surface anisotropy type there appears barrel- or 
pillow- shaped deformation of the vortex core w{z) for 
XT > 0 and xr < 0 respectively."^^ 

Without magnetic held vortices with opposite polar¬ 
ities (p = ±1) are energetically equivalent. Under the 
action of a transversal dc magnetic held the preferable 
polarity of the vortex coinsides with the held direction 
(light vortex). The vortex with opposite polarity (heavy 
vortex) turns into so-called cone phase. Its structure is 
well-known:®®’^^ While the in-plane magnetization dis¬ 
tribution (j) of the cone phase coincides with (6a), the 


out-of-plane component is deformed in the following way. 
The heavy vortex becomes narrower with the typical core 
width Wchih) oc ^J\ — |/i| in contrast to the light vortex 
which becomes broader with Wci{h) oc l/i/l — 

Instead of the exponentially localized structure (6b) 
without held, the vortex prohle has algebraic decay in the 
cone phase. Far from the origin it is characterized by the 
asymptote m = h + hP' For dehniteness we sup¬ 

pose that the vortex polarity p = -1-1 and negative held 
intensity h < 0 having in mind switching phenomenon. 

With the held intensity increasing the heavy vortex 
loses its stability.^® It is important to stress that in spite 
of the fact that the light vortex is energetically preferable, 
the heavy vortex can not perform switching in scope of 
continuous theory due to inhnite barrier which separates 
opposite polarized vortices. However the barrier becomes 
hnite in the discrete spin lattice®® and the reversal pro¬ 
cess can occur®®’®^. The simple picture of the vortex core 
switching provides the core model, originally introduced 
by Wysin®® for the vortex instability phenomenon. Be¬ 
cause of the crucial role of exchange interaction inside the 
vortex core for the switching mechanism, the core model, 
provides good qualitative description of the switching 
process.®"^’®® Let us mention one more possibility, the so 
called cutoff model:®® modeling the discreetness effects 
by cutoff parameter it is possible to describe the vor¬ 
tex polarity switching under the action of spin-polarized 
current®®, by dc^^ and ac^^ magnetic fields. 


III. SPIN-LATTICE SIMULATIONS OF THE 
VORTEX POLARITY SWITCHING 

The main purpose of the current research is to describe 
the fine structure of the vortex core switching process. 
As we discussed above the switching is possible only in 
the discrete lattice. There appear serious difficulties in 
modeling of this process: the standard micromagnetic 
simulators consider the numerically discretized Landau- 
Lifshitz equation, which is valid in continuum theory. 
Such approach becomes mesh-dependent within the con¬ 
tinuum description of micromagnetism.®® The switching 
is known to be accompanied by the creation of micro- 
magnetic singularities, Bloch points, which can not be 
described adequately within the continuum limit. 

In order to overcome difficulties of continuous ap¬ 
proach, we perform spin-lattice modeling using in-house 
developed spin-lattice simulator SLaSi.®® Numerically we 
solve Landau-Lifshitz-Gilbert equations (2) for the disk¬ 
shaped system. We consider two sets of parameters: disk 
diameter 2R = 149a, thickness L = 49a and the mag¬ 
netic length ^ = 14a (disks A), 2R = 99a and £ = 5a 
with the thickness L = {N — l)a, where N is the number 
of lattice sites along disk axis with V = 11, 50 (disks B). 
Chosen diameters allow vortex to be both in the center 
of the sample and between lattice sites. Other parame¬ 
ters are following: the initial vortex polarity p = 1, the 
Gilbert damping constant e = 0.5 corresponds to the 
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overdamped regime. The surface anisotropy is chosen 
x^ES =0.5 and xten = —0.5, and the field intensity varies 
in the range h G (—0.83; 0). The lattice planes are chosen 
parallel to the face surfaces of the samples. 


A. Magnets with bulk anisotropy only: Ks = 0 

Using the vortex distribution ( 6 ) as initial one for dif¬ 
ferent field intensities h < 0, we relax numerically the 
system to the static configuration, which is adopted to 
the lattice. Starting from weak fields, we increase adi- 
abatically its absolute value up to the critical value he- 
For stronger held intensities, \h\ > \hc\, the switching 
occurs. 

We model the switching process for the vortex in 
the disk A. Numerically we found that \hc\ « 0.83. 
The heavy vortex is relaxed in the held with intensity 
h = —0.82. It demonstrates the narrowing of the core 
and shifting of the background far from the origin de¬ 
scribed above. After that we applied the held with in¬ 
tensity h = —0.83 which results in the polarity switch¬ 
ing. The heavy vortex losses its stability: the four cen¬ 
tral magnetic moments in the lattice planes inside the 
vortex core rapidly switch their direction to the opposite 
one. Another magnetic moments relax to the light vortex 
distribution with a smaller speed. In general, the back¬ 
ground of the vortex remains the same, all dynamics is 
observed only in the core region. 

Essentially, the reversal process occurs simultaneously 
in all lattice planes, hence the switching scenario is fully 
two dimensional one. The temporal evolution during the 
switching is shown in the Fig. 2. Schematics (a) and 
(e) represent the 3D magnetization structure in the sam¬ 
ple before and after switching, respectively: vortex state 
within three different horizontal cross-section (top, cen¬ 
ter and bottom surfaces). Central figures (b)-(d) show 
the magnetization distribution within the different cross- 
sections in three different time moments: before the 
switching (r = 6.107), at the switching moment (r = 
6.115), and just after the switching (r = 6.121). During 
the switching process, the linear singularity m{j) = 0 is 
created, see the solid green vertical line 7 in Fig. 2.®^ 

In detailes the polarity switching occurs by the follow¬ 
ing steps, shown in Fig. 3. To consider the spin-lattice 
system we replace the vortex polarity p = ±1 by a dy¬ 
namical polarity /r = /i(z, t) which represents the actual 
amplitude of magnetization ^-component in the vortex 
core center. The magnetic moments in the central part 
of the initial heavy vortex core are approximately perpen¬ 
dicular to the plane and one can say that the dynamical 
polarity is p « 1, see curve (a) in Fig. 3. When the ap¬ 
plied field becomes stronger than the switching threshold, 

rapidly decreases and changes its sign passing through 
fi = 0. One can say about creation of a linear singularity 
in this moment, see curve (b) in Fig. 3. The dynamics of 
the central magnetic moments is going faster than outer 
ones. Central moments rapidly flip their direction from 


/j.R:!lto/r«— 1 , and the overshooting in the vortex 
profile appears, see curve (c) in Fig. 3. The final state is 
the light vortex, see curve (d) in Fig. 3. 

It should be noted that during the whole complicated 
dynamical process of the switching the system remains 
uniform along the z-axis: the numerical difference for 
different lattice planes is not larger than 10 “® for 
(corresponds to the accuracy of the saved data) during 
relaxation in the external field and polarity reversal. The 
same scenario takes place in the case of non adiabatically 
applied external held. 


B. Magnets with the surface anisotropy 


Let us study the role of the surface anisotropy in the 
vortex statics and dynamics, especially in the vortex 
switching phenomenon. Very recently we have studied 
how the surface anisotropy effects on the vortex core 
shape without magnetic held: there appears the pillow- 
and the barrel-deformation of the core for the ES and 
EN anisotropies, respectively.^^ Qualitatively, the vortex 
core width inside the sample volume is determined by the 
magnetic length, while the core width on a surface layer 
is characterized by effective magnetic length"'^ 


4ff 


i 

■\/l -|- 


( 7 ) 


In the case of ES surface anisotropy {k > 0), an effec¬ 
tive magnetic length igg < i. Therefore the vortex core 
becomes more narrow near the face surfaces. This corre¬ 
sponds to the barrel-shaped prohle. In the same way, EN 
surface anisotropy {k < 0) results in the pillow-shaped 
prohle because £eff > see Fig. 4 and schematics (c) 
and (e) therein. 

In weak helds (far from the switching threshold) the 
vortex core looks qualitatively the same: barrel-shaped 
prohle for the ES case and the pillow-shaped one for the 
EN case, while core becomes more narrow in agreement 
with above mentioned estimations. 

Under the action of intermediate held the vortex core 
prohle deforms. The further increase of the held inten¬ 
sity (up to the switching threshold) drastically changes 
the vortex core shape. In the ES case there appears a de¬ 
formed pillow core instead of original barrel one, see the 
Fig. 4 and schematic (f) therein. In the same way in the 
EN case one has the deformed barrel core prohle instead 
of original pillow one, see the Fig. 4 and schematic (d) 
therein. Numerically the vortex core width is computed 
as the wo{z) such that m{wo) = 0.5 [to(0) -I- m{R)]. 

The inhomogeneous shape of the vortex core prohle 
is a key point for understanding the Bloch point medi¬ 
ated switching scenario. Let us start with the case of 
EN surface anisotropy. In this case the vortex core pro¬ 
hle becomes deformed barrel-shaped, see the Fig. 4 with 
schematic (d). With the held increasing both bending 
points shift outside to the ends of the Bloch line and they 
become closer to the face surfaces when held approaches 
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FIG. 2: (Color online) Spatio-temporal picture of the vortex core switching for the magnet with bulk anisotropy 
only, the simulation data (center). Left (a) and right (e) panels: schematic of the vortex profiles in three different 
horizontal cross-sections. During the switching process the vortex polarity flips simultaneously in all lattice planes. 
The gray rectangle on background indicates the XZ cross-section of the sample and the white boxes show the spatial 
profile of the ^-component of magnetization mz(r, z). The magnetization distribution is shown in three different 
lattice cross-sections (zt = 24.5a, Zc = 0.5a and Zh = —24.5a) at three different time moments: b) before the 
switching, r = 6.107, c) at the switching moment, r = 6.115, and d) just after the switching, r = 6.121. The 
switching is accompanied by the linear singularity creation, see the green line. Parameters: disk A, h = —0.83. 



Distance r/a 


FIG. 3: (Color online) Vortex polarity reversal in the 
sequential time moments: (a) Initial state of the heavy 
vortex before switching; (b) Dynamics polarity in the 
center decreases due to discreetness, linear singularity 
appears; (c) Central magnetic moments reverces faster 
than other ones, overshooting appears; (d) Final state 
of the light vortex. 


to the switching one. This means that in a strong enough 
field the vortex profile near the face surface is more nar¬ 
row than in the bulk. 

This bulk deformation of the vortex core results in the 
inhomogeneous vortex polarity switching: the Bloch line 
breaks in its more narrow parts, i. e. on the face sur¬ 
faces. Each break of the Bloch line corresponds to the 
Bloch point. Therefore in case of EN surface anisotropy 
two Bloch points are nucleated simultaneously on the face 
surfaces. During the reversal process they move in op¬ 


posite direction to each other and finally annihilate at 
the sample’s center. The similar behaviour occurs in the 
case of ES surface anisotropy: the narrowest part of the 
vortex core is situated in the center and a pair of Bloch 
points is nucleated in the center of the sample. 

The schematics of the 3D polarity reversal accompa¬ 
nied by Bloch points motion for EN surface anisotropy is 
shown in Fig. 5. The simulation data are shown for the 
disk A under the action of the external field h = —0.85 
applied to the sample relaxed in h = —0.8. Bloch points 
always separate vortices with opposite polarities. There¬ 
fore their existence can be easily detected by observing 
of the vortex polarities on the different lattice planes: 
face surfaces and in the center of the sample, see Fig. 5a. 
The temporal evolution of both Bloch points has the re¬ 
flection symmetry with respect to z = 0 plane: curves 
which show polarity dynamics of top and bottom vor¬ 
tices coincide (solid curve). The switching of the central 
vortex occurs with delay (dashed curve). In thicker sam¬ 
ples more pairs of Bloch points can be nucleated inside 
the sample volume, see below (also see inset in Fig. 7). 

The axial cut of the lattice during polarity switching is 
shown in the Fig. 5b: the vortex state is already reversed 
(/i « —1) near the top and bottom surfaces, while its po¬ 
larity is still directed up (/r « -1-1) in the sample’s center. 
There are two Bloch points in the system: their positions 
are shown by thick white arrows and the direction of mo¬ 
tion is shown by thin green arrows. The vortex prohles 
for different times during polarity switching process are 
shown in the Fig. 5c-e. 

Let us discuss now the case of ES surface anisotropy, 
see Eig. 6 (notations are the same as in Fig. 5). Here we 
apply the external field h = —0.85 to the sample relaxed 
in h = —0.8. It results in the deformed pillow-shaped 
vortex core profile. That is why the vortex core becomes 
more narrow in the middle of the sample, hence the Bloch 
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Dimensionless vortex width w 



Dimensionless vortex width w 


(a) EN surface anisotropy, = —0.5 


(b) ES surface anisotropy, = 0.5 


FIG. 4: (Color online) The reduced vortex core width ic as a function of the sample thickness coordinate z for the 
disk B with thickness L = 39a. Without field the vortex core has the pillow shape for the EN case (xr = —0.5), see 
schematic (c), and the barrel one for the ES case (xr = 0.5), see schematic (e). Under the action of the field close to 
switching one (h = —0.6) there is an opposite behavior: the vortex core profile becomes a deformed barrel one for 
the EN case, see schematic (d), and a deformed pillow one for the ES case, see schematic (f). Arrows show 
transformation under the action of external field h which is increased by modulus. 



EIG. 5: (Golor online) Switching of the vortex polarity for 3D Heisenberg magnets with EN surface anisotropy 
(disk A, x: = —0.5, simulations data), a) Polarity dynamics for vortices with z = 24.5a and z = —24.5a (solid curves 
which coincide), z = 0.5a (dashed curve), b) Magnetization distribution in the x-z plane near the center of the 
sample at r = 9.07. Color indicates the m sign, thick arrows show the positions of the Bloch points, thin vertical 
arrows indicate the direction of motion for the bottom and top Bloch points, c-e) Polarity switching in different x-y 
planes, according to curves on a). Green arrows show the direction of the Bloch points (BP) motion. 


line breaks in the middle and the switching process starts 
inside the sample. There are two Bloch points which born 
in the sample center. During the switching process they 
move along z-axis in opposite direction to face surfaces. 

Numerically we also analyzed the case of the exchange 
anisotropy; qualitatively the mechanism is very similar 
to the EN surface anisotropy, see Appendix B for details. 


To gain insight into the temporal evolution of the 
switching process we performed a set of simulations for 
fixed diameter of the sample (disk B) and varying thick¬ 
nesses {N = 11,40). The anistropy is xr = —0.5. We 
apply the field h = —0.8®® to the vortex state relaxed 
without field {h — 0). Similarly to the previous case, the 
inhomogeneous polarity reversal is observed. In all sim- 

























































































FIG. 6: (Color online) Switching of the vortex polarity for 3D Heisenberg magnets with additional ES surface 
anisotropy {x = 0.5), the simulations data. All notations and other parameters are the same as in Fig. 5. 




(b) ES surface anisotropy, x: = 0.5 


FIG. 7: (Color online) Bloch point position as function of time for samples with different thickness. All curves are 
matched at r = 0 for EN surface anisotropy and at r = 0.548 for ES surface anisotropy. Green circles indicate Bloch 
points creation event and crosses indicate annihilation event. Eor N > 30 additional pair of Bloch points is 
nucleated shortly before annihilation. Arrows indicate direction of motion of Bloch points. Parameters: disk B, 

X = —0.5, h = —0.8. 


ulations the reversal is accompanied by the nucleation of 
a pair of Bloch points near face surfaces, their further 
dynamics inside the magnet and the final annihilation at 
the point of contact in the middle. The temporal dynam¬ 
ics of the Bloch points for different thicknesses is shown 
in Eig. 7. We determine the position of Bloch points 
as a cross-sections of three isosurfaces ma{r) = 0 with 
a = x,y, Each surface ma(r) is computed as interpo¬ 
lated functions for magnetization components based on 
discrete values of m„. 

A new feature appears for relatively thick samples 
{N > 30): During the switching process, additional pair 
of Bloch points is nucleated in the sample center. Two 
points of this pair repel each other moving in z-direction 
out from the center. Einally the new born Bloch points 
annihilate with originally nucleated ones, see the inset in 
Fig. 7. Note that for the weaker field h = —0.7 the addi¬ 


tional Bloch points are not nucleated. One can conclude 
that there appears instability of the Bloch line in long 
enough samples, which causes the breaking of the Bloch 
line inside the sample and nucleation of additional Bloch 
points. Such a picture is similar to the temporal forma¬ 
tion of micromagnetic drops in long nanocylinders. 


Numerically we also studied how the speed of the Bloch 
point depends on the sample thickness. One can see from 
the Fig. 7 that the Bloch point motion is almost a steady- 
state one except the moments of nucleation and annihila¬ 
tion. In this way we computed the speed of Bloch points 
in a wide range of the sample thickness, see Fig. 8. 
The Bloch point speed rapidly decreases in thick sam¬ 
ples and tends to some constant value not dependend on 
the sample thickness and surface anisotropy type. 
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FIG. 8: (Color online) Bloch point speed as function of 
sample thickness in units of lattice sites along disk axis 
for the EN and ES surface anisotropies. Parameters are 
the same as in Fig. 7. 


IV. DESCRIPTION OF THE POLARITY 
SWITCHING PROCESS 

To describe analytically the observed vortex polarity 
switching process we propose a wired core model. In 
spite of simplicity this model is in a good qualitative 
agreement with full-scale simulations. It generalizes the 
discrete reduced vortex core model initially proposed by 
Wysin for the vortex instability phenomenon and later 
extended for the description of the polarity switching 
in the two-dimensional case.^®’®°’®^’®®’®^ In the current 
study we extend the model for 3D systems. A wired 
core model allows one to take into account processes of 
nucleation-annihilation of point singularities during the 
vortex switching. 

We start with the discrete Hamiltonian (1). The 
anisotropy term (Ib) can be written in the following form: 

r " _ ( 8 ) 

^ \K, forn = 2,V-l 

" + Ks, for n = 1 and n = N. 

Here n = {nx^riy^n) and n enumerates the planes along 
z direction. 

We suppose that in each n-th plane there exist only 
four “free” magnetic moments, all other moments are 
“frozen” into the background, see Fig. 9. The 3D vortex 
core is considered in our model as a wire of free magnetic 
moments. 

In original core model^® the discrete core magnetiza¬ 
tion is matched with the continuous magnetization of the 
background (the ring), where the magnetization is fixed 
in the pure planar vortex configuration. The inner ra¬ 
dius Tin of the ring should be chosen from the condition 


that the ground state is a pure planar vortex, which ful¬ 
fills when Tin > 0.3£.®® The situation drastically changes 
under the field action, the planar vortex can never be 
a solution of the ring wire in the field presence. Besides 
one has to take into account z-dependence of the magne¬ 
tization due to the surface anisotropy in the way similar 
to (5). Finally, we get the background magnetization 
distribution TOr(r, z), see Appendix C for details: 


= mb{Zn) = hmr{r2,Zn), 

TT (9a) 

4’{n^,ny,n) “ X + ) |^a;| + l^yl > 2, 

where r 2 is the radius of the second coordination sphere, 
function mr is defined by Eq. (C3), z„ is a z-coordinate 
of n-th lattice plane. For the numerical investigation of 
the wired core model we use r 2 ~ 0.35f. 

All “free” moments come from the ring of principal 
sites, they form the vortex core. Due to the axial sym¬ 
metry, all free moments have the same out-of-plane com¬ 
ponent and the same in-plane phase tp which describes 
the deviation from the vortex distribution (9a): 

'^(n^,ny,n) ~ ^ 

TT ( 1 ^^) 

4^{nx^ny.n) X d” ~ -f '0n(O; 'Ox: 'Oy 1, 1. 

We consider the dynamical vortex polarity /iri(t) for each 
lattice plane, and the in-plane magnetization angle 'ipn(t), 
which has the meaning of the turning phase, as two col¬ 
lective variables. In this manner we generalize the 2D 
approach, which we recently proposed in Ref. 15. Note 
that the dynamical vortex polarity is mutually connected 
to the effective vortex width in the full-scale model. 

Now by incorporating the wired core Ansatz (9) into 
the Hamiltonian (1) with account of the anisotropy term 
(8), one can write down the model Hamiltonian, normal¬ 
ized by K§^ as follows 


N 

^c=4^ 

n—1 

4A 


(1 -|- — 2A) — {h 2Am{,) Hr, 


75 


{1 - ml){l - p.1) cos 


N-l 


n—1 


X COs{l{jn - tpn+l) ■ 


( 10 ) 


■ 4A [l^nlJ-n+l + \/(l-M^)(l-/i^+l) 


Here = 0 for n = 2,N —1 and = >«■ for n = 
1,N, see Appendix D for details. The first sum in (10) 
describes the intra-plane interaction and the second sum 
describes the inter-plane interaction. 

The equilibrium values of Pn and ipn can be found 
analytically for the homogeneous case Xn = 0, n = 1, N 
and h — 0: 


phom{h — 0) — it 



16A2 

5(2A- 1)^’ 


ipo{h = 0 ) = 0 , 
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Bloch line 


FIG. 9: (Color online) Schematic of the wired core 
model for n-th and (n + l)-th planes: Red arrows 
indicate free magnetic moments (Eq. (9b)) and gray 
ones indicate fixed ones (Eq. (9a)). Layers are 
exchange-coupled with exchange integral J and the each 
plane is characterized by own anisotropy constant 
Kn > 0, see Appendix D for details. Bloch point 
position is marked by the green sphere. 


where A > 10 for the out-of-plane vortex. For h ^ 
the dependence /ihom(h) can be found numerically as the 
solution of the Eq. (D3), see Eig. 15 in Appendix D. 

When the field is applied one has to take into account 
that the equilibrium z-components of magnetization are 
different for the bulk and surface planes. Eor the hxed 
magnetic moments the out-of-plane magnetization com¬ 
ponent is taken equal to mb{z) according to Eq. (C3). 
Assuming that the deviation x{z)>clV\ = fi{z) — /ihom 
from the equilibrium value fihom is small, <C 1, 

and replacing the Hamiltonian (10) by its continuum ver¬ 
sion we obtain the following boundary value problem for 
x{z): 


a x”{z) - Aix{z) = A 2 (z), 


±ax'(z) + 


Mhom (1 Mhom ) 


Ta 


J=±L/2 


= 0 , 


( 11 ) 


where coefficients Ai and ^ 2 ( 2 ) are given by Eqns. (D4), 
see Appendix D for details. Function A 2 (z) contains only 
simple harmonics and hyperbolic sine and cosine, so the 
solution of Eq. (11) can be found analytically, neverthe¬ 
less it is unwieldy and below we present only numerical 
solution. 

The dynamical polarity as numerical solution of 
Eq. (11) for both cases of surface anistoropy is plotted in 
Eig. 10 by solid lines. One can see this result is in a good 
agreement with a direct solution for Hamiltonian (10) 
(symbols correspond to simulations of the system with 


ES, model • EN, model ■ 



fJ^n(z) 


EIG. 10: (Color online) Dynamical polarity /i„ for the 
discrete model (10) (symbols labeled EN and ES model) 
and fj,(z) as numerical solution of Eq. (11) (solid curves 
labeled EN and ES linear approx.). Parameters: 

|>f| = 0.1, A = 20, A^ = 70, h = -0.0009. Dashed line 
corresponds to polarity without surface anisotropy. 


|>!:| = 0.1, A = 20 and N = 70). Numerically we min¬ 
imize (10) with respect to /x„ and ipn under the action 
of held h = —0.0009. In equilibrium tpn = 0. The cor¬ 
responding equilibrium polarity in the homogeneous sys¬ 
tem /ihom(h) = 0.266 is shown by the dashed line. 

The dynamical polarity in the wired core model plays a 
role of the vortex width w{z) in the full-scale simulations. 
They reproduce the same shapes near the switching held, 
cf. Eig. 4 and 10, and the switching starts in places with 
minimal /i„. The deformation of the vortex shape in the 
wired core model is a result of the inhomogeneous ground 
state mb{z) and is highly sensitive to its exact form. 

The switching occurs when held increases to \h\ = 
0.001. As in the full-scale simulations discussed above 
the Bloch points nucleate on the surface for EN surface 
anisotropy, see Fig. 11a, cf. Fig. 7a; they propagate from 
the lattice planes with n = 1 and n = N to the center of 
the axis. For a thick enough lattice systems additional 
Bloch points are nucleated, see curve for N = 140 in the 
Fig. 11a. Bloch points are nucleated in the center of the 
system for ES surface anisotropy in the same way as in 
the full-scale simulations, see Fig. 11b, cf. Fig. 7b. 

We also studied thickness dependence of the Bloch 
point speed in the same way as in full-scale simulations, 
see the Fig. 11 for the held h = —0.001 (numerical so¬ 
lution of Eqns (Dl) with e = 0.5. The comparison with 
the full-scale simulations (cf. Eig. 8) conhrm the fast de¬ 
cay of the Bloch wall speed with increasing of the sample 
thickness with the further saturation. 
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(a) EN surface anisotropy, = —0.1 



(b) ES surface anisotropy, x: = 0.1 


FIG. 11: (Color online) Bloch point position as function of time for samples with different thickness for the wired 
core model. All curves are matched at t = 0 for EN surface anisotropy and at t = 0.125 for ES surface anisotropy. 
Green circles indicate Bloch points creation event and crosses indicate annihilation event. For N = 140 additional 
pair of Bloch points is nucleated shortly before annihilation. Arrows indicate direction of motion of Bloch points. 
£ = 0.5, h = —0.001, other parameters are the same as in Fig. 10. 



Thickness N 

FIG. 12: (Color online) Bloch point speed as function of 
sample thickness in the wired core model, cf. Fig. 8. 
Parameters are the same as in Fig. 11. 


V. DISCUSSION 

To summarize, we described different mechanisms of 
the vortex polarity switching in Heisenberg magnets with 
and without additional surface anisotropy. Under the 
action of transversal dc magnetic field the vortex core 
is reversed using the axially-symmetric scenario. 

We conclude that without surface anisotropy the switch¬ 
ing occurs uniformly with respect to the thickness z- 
coordinate through the transient linear singularity or the 
planar vortex. In presence of the surface anisotropy the 
switching is accompanied by nucleation of point singular¬ 
ities (Bloch points), their motion and final annihilation. 


The birthplace of Bloch points and the direction of their 
motion depends on the type of the surface anisotropy. 
The complicated vortex dynamics including its switch¬ 
ing obtained using full-scale spin-lattice simulations can 
be described analytically by a wired core model, which is 
elaborated in this work. 

Qualitatively, the influence of the surface anisotropy 
can be explained as follows. In terms of the surface 
anisotropy the magnetization configuration is pinned 
near the sample boundary. This effective pinning has to 
compete with exchange interaction, which results in the 
Robin type boundary condition (4b). That is why the 
homogeneous magnetization distribution m = h is not 
possible, since it does not satisfy the boundary condi¬ 
tion; there appears mo(z,h) background profile, see (5). 
Similar scenario takes place for the vortex state particle, 
where the vortex width becomes ^-dependent, w = w(z). 
According to (4b) the sign of the vortex core width gra¬ 
dient depends on the surface anisotropy type: barrel¬ 
shaped for the ES surface anisotropy and pillow-shaped 
for the EN one.^^ 

Let us sketch the physical picture of the influence of 
the surface anisotropy on statics and dynamics using the 
particular case of ES surface anisotropy. The typical vor¬ 
tex core width ic on a surface layer without magnetic 
field is determined by effective magnetic length £eg, see 
(7). In ES case, the value of effective magnetic length 
£eff < £i and the vortex core becomes more narrow near 
the face surfaces, which corresponds to the barrel-shaped 
profile. Under the action of magnetic field the vortex 
state undergoes several changes. Eirst of all, there ap¬ 
pears a cone state instead of the easy-plane one. Apart 
this the core width for the heavy vortex becomes nar¬ 
rower, wiji) cx \J\ — \h\. However there is a counter¬ 
action of the surface anisotropy, which tries to fix the 
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(a) EN surface anisotropy, nanodisk of radius 100 nm 
and thickness 24 nm 



\h\ grows 




(b) Effective ES surface anisotropy induced by dipolar 
interaction, nanodisk of radius 66 nm and thickness 
100 nm 


FIG. 13: (Color online) Simulations with dipolar interaction using OOMMF for EN and effective ES surface 
anisotropies. Schematics show the change of the vortex prohle in the same way as in Fig. 4 


z-gradient of the core width w{z). As a result of this com¬ 
petition, there appears a neck in the vortex core profile, 
which results in the deformed barrel core instead of orig¬ 
inal barrel one, Fig. 4 with schematics (e) and (f). This 
static picture plays a crucial role for the understanding 
the reversal mechanism. When the held intensity ap¬ 
proaches the critical value (switching held), the vortex 
core width extremely decreases. However this happens 
nonhomogeneously: the Bloch line has a bottleneck in 
its center, hence it breaks during the switching in this 
place. Namely this point is the birthplace of two Bloch 
points, which move along the Bloch line repelling each 
other until they annihilate on the face surfaces of the 
sample. 

In the current study we do not take into account the 
dipolar interaction, which is of great importance for the 
real magnets. Specifically, namely the dipolar interac¬ 
tion favors the magnetization curling and causes that 
the vortex becomes the ground state of magnetically soft 
nanodisk.In our previous work on the influence of the 
surface anisotropy on the vortex core we compared in de¬ 
tails different limiting regimes when magnetostatics can 
lead to effective anisotropy, see Ref. 44 and references 
therein. 

In order to validate our results for nanomagnets with 
the dipolar interaction, we also performed few simula¬ 
tions with account of both surface anisotropy and dipo¬ 
lar interaction. In the current study we do not use SLaSi 
simulator for the dipolar interaction case, because of 
high numerical costs. Instead we make micromagnetic 
OOMMF^^ simulations. As we already discussed above 


the modeling of the Bloch point is challenging due to 
singular magnetization distribution, leading typically to 
the mesh size dependence.Nevertheless, according to 
the physical picture of the described effect, the key role 
in the dynamic switching phenomenon plays the static 
instability of the vortex core: there appears a spatial 
nonhomogeneity of the Bloch line, which causes breaks 
of Bloch line at the bottleneck places and finally results 
in the Bloch point mediated switching. That is why we 
perform OOMMF simulations of different static vortex 
distributions which precede switching. Numerically we 
model the Permalloy sample with account of dipolar in¬ 
teraction without surface anisotropy (because the magne¬ 
tostatics produces the effective ES one"’’^) and additional 
EN surface anisotropy.®® 

Similar to the Heisenberg magnet with EN surface 
anisotropy the vortex core in the nanodisk with magne¬ 
tostatics and strong enough EN surface anisotropy has a 
pillow-shaped profile without external field, see Fig. 13a 
[data for FI = 0 and schematic (c)]. Under the action of 
the field there appears the curve bend of the vortex core 
profile. Under the action of strong enough field the vor¬ 
tex profile takes up a deformed pillow shape and barrel 
shape in the stronger field just before switching field, see 
Fig. 13a [data for H = —3 kG (0.3 T) with schematic (d) 
and data for H = —4.4 kG (0.44 T)]. Even in the absence 
of additional ES surface anisotropy in the sample with 
large aspect ratio L/{2R) > 1 the initial barrel-shaped 
vortex core deforms in the central part in the same case 
as it is shown in Fig. 4e. Such a bending is absent in the 
thin samples where the creation of volume magnetostatic 
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charges significantly infiunces to the surface ones. In the 
both cases of surface anisotropies switching starts from 
the places, where vortex core becomes more narrow un¬ 
der the action of external field: from the surfaces for EN 
and from the center for effective ES surface anisotropy. 
Thus, the influence of the external field on the vortex 
profile and switching mechanism for the EN and ES sur¬ 
face anisotropies is in a good qualitative agrement with 
results of SLaSi simulations, see Fig. 4 and wired core 
model, see Fig. 10. 
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Eq. (B) runs over 6 neighbors for volume sites and over 3- 
5 for surface sites which can be interpreted as a exchange- 
induced surface anisotropy due symmetry breaking of the 
lattice on the surface. 

We use the same procedure of simulations as in previ¬ 
ous sections with JzjJ = 0.005, which gives the vortex 
polarity reversal under the action of the reduced mag¬ 
netic field h = 2 hbH/{J z§) = 2.2. The observed magne¬ 
tization dynamics is qualitatively similar to the case of 
easy-plane volume anisotropy with EN surface anisotropy 
shown in Fig. 5. In comparison with simulations with 
single-ion anisotropy the difference in vortex core mag¬ 
netization is much smaller in this case. We compute the 
difference of dynamical polarities in the center and face 
surfaces of the sample = |/itop ~ Mcenterl = 1-65 

(see Fig. 5d) and = 0.032 which is about 50 times 

weaker than EN surface anisotropy. 


Appendix A: Equilibrium magnetization in film with 
surface anisotropy 

Let us consider a film of thickness L with constant of 
the easy-plane anisotropy K > 0 and constant of the 
surface anisotropy = kK. The equilibrium magne¬ 
tization distribution is determined by the Eqns. (4). In 
order to derive the magnetization distribution one can 
assume that m = \/\ — vri^x -|- moz without lost of gen¬ 
erality, where mo = mo{z,h). In this case Eqs. (4) are 
reduced to 


Appendix C: Ring under the action of DC magnetic 

field 

Let us consider a magnetic ring of inner radius > 
0.3f in the vortex state described by the Hamiltonian (1). 
For such an inner radius without external fields the out- 
of-plane magnetization component vanishes and planar 
vortex appears.®® The Eq. {4a) takes the following form 
after substitution of the vortex anzats (6a): 

2 

771 ^ 

i'^V^niz + h — ruz = —^V^m„, (Cl) 

niy 


2 d'^mo 
dz^ 


+ h — mo 


£^mo 9^ 
Vl - mg 



'fdmo 

dz 





= 0 . 

z=±L/2 


(Al) 

(A2) 


The analytical solution can be found in assumption that 
\x/'</\\ 1. Assuming mo(z, h) = h + -^y(z) the solu¬ 

tion takes the form (5). 


Appendix B: Simulations with exchange anisotropy 


The surface anisotropy naturally appears in the 
Heisenberg ferromagnet with exchange anisotropy. 
That is why we additionally simulate the system, de¬ 
scribed by the LLG equations (2) and Hamiltonian 


n,d 

- 2fiBH§ E 




E 

n.S 


TR-nUlri-\-S 


where niy = y/l — m^ cos (j). Let us consider the case of 
weak fields \h\ <C 1. We find the solution of Eq. (Cl) in 
the form 


rUz = hrrir. 

The Eq. (Cl) takes the form 
' d^nrir d^nir ^ 1 dm 


dz^ dr"^ r dr 


dmr xa 


dmr 


dr 


Wr I ^ - 1 


= 0 , 


z=±Ll2 


= 0 . 


-bl = 0 
(C2a) 

(C2b) 


(C2c) 


r=rin,i? 

For >r 0 the solution of (C2) can be found in the form 
m^ = ^ Rn{r)Zn{z) (C3) 


which results in the following boundary-value problem: 


where J 2 > 0 is a coefficient of the exchange anisotropy 
for easy-plane magnet. The sum in the second term in the 



+ — 0 , 


(C4a) 
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FIG. 14: (Color online) Magnetization distribution 
m,z = hmr{r, z) in the center and on the face surfaces of 
the ring with rin = 1.58a and ES surface anisotropy 
under the action of the field h = —0.0009 (4 harmonics 
of the Eq. (C4) are taken into account). Dotted curve 
corresponds to the function hnirh (absence of surface 
anisotropy). Inset shows value of coefficient B as 
function of po = \/\ + = 0.1, other 

parameters are the same as in Fig. 10. 


d^Rn 1 di?„ 


dr^ r dr 


Rn I ~ ^ ~ ] d- Fn — 0^ 


dz^ 

dz ^ Vx ^ 


= 0 , 


z=±L/2 


(C4b) 

(C4c) 


di?„ 


dr 


= 0 . 


(C4d) 


r—Tin, ft 

where Fn are defined as coefficients of expansion 

l = ^F„Z„(z). (C5) 


Eigenfunctions Zn have the form 

crn{z + L/2) XT , (7„(z + L/2) 

Zn{z) = cr„ cos ^ sm -- J 


with dispersion equation 
^ (7 ri.L 




For the case >r < 0, the first eigenvalue cto should be 
found as solution of 


tanh 


ctqL 

2£ 


VXao 


with eigenfunction 


ry / ^ ^(^0{z + L/2) XT , (7o{z + L/2) 

Zo{z) = (To cosh- j -h ^ smh- - -. 

Solution of the Eq. (C4b) can be written in the following 
form 


Rn{p) =-- 


Fr,. 


Kj(/o) / L^{p)pdp + BKi{p) 


OO 

+hi[x) J Ki{p)pdp 


(C6) 


where p = ^1 + (T^r/£, Li(*) and Ki(*) are real val¬ 
ued solutions of modified Bessel equation of purely 
imarginal order with Wronskian 'W{Ki((C),L(C)} = 

2 /c 73,74 Constant B can be found from the boundary 
conditions (C4d) at r = 


B{po) 


U{p)pdp- 


ddjj/dp 

dKi/dp , 


Ki{p)pdp 


B^\po) 


Po +0.557 
0.99po - 0.316’ 


P=P 0 


(C7) 


where po = a/I + The interpolation function 

B^^{po) fits numerically calculated dependence B{po) in 
the range po € [0.34,0.8] with an accuracy of about 
5%. Functions mr(r) and B{po) are shown in Fig. (14). 
In the absence of the surface anisotropy the solution of 
Eqns. (C2) can be found in the form rur = mrhii") which 
is the solution of the following equation: 



I dmrh 

r dr 


TTlrh 


S-1)+1 = 0 (C8) 


with boundary condition (C2c). Its solution can be ob¬ 
tained from Eq. (C6) by replacing to nirh and sub¬ 
stitution Fn = 1 and (t„ = 0 . 

For the calculations in scope of the reduced wired core 
model we use 3 harmonics for smallest L and up to 7 for 
the most thick samples. 


Appendix D: Reduced Wired Core Model 

The Hamiltonian (1) with account of the anisotropy 
term (8) can be rewritten in the following form: 

Htf, — JS ^ ^ 1^777^777^^5 
n,S 

+ cos((/)„ - (l>r,+s) 

+ . 

n—{nx ,ny ,n) 
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Field 10^ h 


their continuous analogues fj,{z) and ipiz) and replace 
summation by integration. The volume contribution has 
the form 


L/2 


(1 — 2A) —— {h + 2Amb)/i 


-L/2 


4A 

■y! 



L/2 

-|- 2Aci / 


-L/2 


/\2 




dz 




(D2a) 


dz, 


FIG. 15: (Color online) Equilibrium polarity /rhom(h) 
for the discrete homogeneous model (10) (xr^ = 0, 
n = I, N) for A = 20. 


By incorporating the reduced vortex core ansatz (9) one 
get 


N 
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- 4J§2 + Y^(l-M^)(l-Mn+l) 


n=l 


with Too = wo(z, h) defined in Eq. (5). The surface con¬ 
tribution is 




2xr/r^ 


z=±Ll2 


(D2b) 


We assume that the deviation from the equilibrium dy¬ 
namical polarity in absence of the surface anisotropy is 
small and ^(z) = ^o(^) + :^ 2 :(z) with \x/'/X\ 1, 

where fJ,hom{h) can be found numerically as solution of 
the following equation 


2 (1 2A)^j^qj^ /l(l -f 2A77Z7'/i)/iho 

4A 




(D3) 


X cosiljjn - 1pn+l) 


Taking into account that >«■„ = 0 for n = 2,iV — 1, >!:i = 
kn = XT, a = P ja^., A„ = JjKn = A/(l -I- x:„) one 
obtains the Eq. (10). The temporal evolution of ^nij) 
and '0n(T) is governed by the Eq (2), which results in 


see Fig. 15. Taking ip{z) = 0 and varying the func¬ 
tional (D2) one can find the equilibrium distribution of 
/i as solution of the boundary-value problem (11) where 

^1 = ( 2 -(l-/r^om) -(D4a) 

V V5(l-FLm) 


du„ 1 dM’c 2\ £ dJfc 

dV>„_ I dM’c 1 £ dJ^c ^ ^ ^2(2:) 

dr 1+ d/j,n 1 - 1 -|- ' 

The problem of magnetization dynamics in discrete 
model can be analysed numerically. For analytical ap¬ 
proach we replace the discrete functions Un and ijir, using _ 
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+ Mh 


om 


1^5 - lOA -f 4A 



(D4b) 
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